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Abstract 

Q-i We present exact solutions to the Einstein-Maxwell system of equations in 

5^ ' spherically symmetric gravitational fields with a specified form of the electric 

field intensity. The condition of pressure isotropy yields a difference equation 
with variable, rational coefficients. In an earlier treatment this condition was 
integrated by first transforming it to a hypergeometric equation. We demonstrate 



> 

in 

CN ■ that it is possible to obtain a more general class of solutions to the Einstein 

£T) • Maxwell system both in the form of special functions and elementary functions. 

Our results contain particular solutions found previously including models of 
^ ■ charged relativistic spheres and uncharged neutron star models. 

o 



1 Introduction 



Exact solutions of the Einstein-Maxwell system are important in the description of rel- 
ativistic astrophysical processes. In the presence of charge the gravitational collapse of 
a spherically symmetric matter configuration to a point singularity may be avoided [l|. 
Einstein-Maxwell solutions are important in studies involving the cosmic censorship 
hypothesis and the formation of naked singularities 0]. The presence of the electro- 
magnetic field affects the values of the surface redshifts, luminosities and the maximum 
mass for charged stars as demonstrated by Sharma et al [3|, Ivanov 4| and others. The 
analysis of Patel and Koppar [5], Patel et al [sj, Tikekar and Singh 7|, Gupta and Ku- 
mar jg], and Mukherjee 9| have shown that it is feasible to model charged superdense 
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neutron stars with densities in the range of 10 14 g cm -3 ; these models generate bounds 
for the surface readshift, luminosity and stellar mass which are consistent with ob- 
servation. Sharma et al [lCj studied cold compact spheres, Tikekar and Jotania 11] 
analysed strange matter, Sharma and Mukherjee [3] considered the equation of state 
of a compact X-ray binary pulsar Her X-l, and Sharma and Mukherjee jl3| investigated 
stars composed of quark-diquark particles which are consistent with the presence of an 
electromagnetic field. Einstein-Maxwell solutions may be applied to the core envelope 
models of Thomas et al [jjj], Tikekar and Thomas 15], and Paul and Tikekar 16] 



with an inner core surrounded by an outer layer. These references provide a sample 
as to why the Einstein-Maxwell system, describing the interior of a charged star, has 
attracted the attention of many researchers. 

Our intention in this paper is two-fold. Firstly, we seek to model a charged rel- 
ativistic sphere which is physically acceptable. We require that the gravitational, 
electromagnetic and matter variables are finite, continuous and well behaved in the 
stellar interior, the interior metric should match smoothly with the exterior Reissner- 
Nordstrom metric, the speed of sound is less that the speed of light, and the solution is 
stable with respect to radial perturbations. Secondly, we seek to regain an uncharged 
solution of Einstein equations which satisfy the relevant physical criteria when the 
electric field vanishes. This ensures that a neutral relativistic star is regainable as a 
stable equilibrium state. We seek to construct a model, with limiting uncharged stars 
as known exact solutions, that exactly satisfies the Einstein-Maxwell systems. This 



ideal is not easy to achieve in practice and only a few examples wit 
features have been found thus far. Recently Hansrai and Maharaj 



i the required two 



171 ]. Thirukkanesh 



and Maharaj [181 ]. and Komathiraj and Maharaj |l9l. |20| found charged relativistic fluid 
spheres which regain neutral compact stars for particular parameter values. 

Our objective is to find new solutions of the Einstein-Maxwell system that satisfy 
the physical criteria given above, and necessarily contains a neutral stellar solution. 
Our approach here complements the approach of Thirukkanesh and Maharaj [l8| who 
presented exact solutions to the Einstein-Maxwell equations; this family of solutions 
can be written in terms of elementary functions and reduces to the well known un- 
charged stellar model of Durgapal and Bannerji 21|. The approach in 18| introduces 
a transformation that reduces the condition of pressure isotropy to a hypergeomet- 
ric equation. The transformation utilised in 18j does produce new exact models but 
restricts the classes of solutions that are possible because of constraints placed on 
particular parameters. In this paper we do not transform the condition of pressure 
isotropy to the hypergeometric equation but are still in a position to integrate the field 
equations. A new class of Einstein-Maxwell solutions are found that contain familiar 
uncharged models which are regainable for different choices of the metric function and 
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the electric field. In Section 2, we express the Einstein-Maxwell field equations for 
the static spherically symmetric line element as an equivalent set of differential equa- 
tions utilising a transformation from 21]. We choose particular forms for one of the 
gravitational potentials and the electric field intensity. This enables us to obtain the 
condition of pressure isotropy in the remaining gravitational potential which facilitates 
the integration process. In Section 3, we assume a solution in a series form which 
yields recurrence relations, which we manage to solve from first principles. It is then 
possible to exhibit exact solutions to the Einstein-Maxwell system. Solutions in terms 
of elementary functions are possible for particular parameter values. In Section 4, we 
present two linearly independent classes of solutions as combination of polynomials 
and algebraic functions. In addition we show that it is possible to express the general 
solution of the Einstein-Maxwell system in terms of elementary functions. We regain 
some solutions found previously from our general solutions in Section 5. In Section 
6, we summarise the results found and present plots illustrating the behaviour of the 
matter variables. 



2 The isotropic equations 

We take the line element for static spherically symmetric spacetimes to be 

ds > = - e Mr) df 2 + e 2Mr) dr 2 + ^2 + ^2 (1) 

where u(r) and A(r) are arbitrary functions. For charged perfect fluids, the Einstein- 

P+\E 2 (2a) 
V ~ \E 2 (2b) 
P+\e 2 (2c) 

y 2 e~\r 2 E)' (2d) 

for the geometry described by ([1]). The energy density p and the pressure p are mea- 
sured relative to the comoving fluid 4- velocity u a = c~ u 5q, E is the electric field in- 
tensity, a is the proper charge density, and primes denote differentiation with respect 
to r. In the system (I2al) - fl2d|) we are using units where the coupling constant = 1 
and the speed of light c = 1. The field equations (!2a|) - pd|) govern the behaviour of 



Maxwell field equations can be expressed as follows 

i[r(l-e- 2 *)]' = 

4(1-0 + — e^ 2A = 

-2A f H , 12 V 1 , , X'\ 

e \ v + v H v A = 

\ r r J 
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the gravitational field for a charged perfect fluid. When E = we regain Einstein 
equations for a neutral fluid. 

An equivalent form of the field equations is generated if we introduce the transfor- 
mation 

A 2 y 2 (x) = e 2u{r) , Z(x) = e~ 2A(r) , x = Cr 2 (3) 

where A and C are arbitrary constants. Under the transformation Q, the system ([2]) 
becomes 

(4a) 
(4b) 
(4c) 





1 - Z 

2Z = 

X 


p E 2 
C + 2C 




y Z - 1 
4Z y + 

y x 


p E 2 
C~2C 


AZx 2 y + 2Zx 2 y + 


( ■ E 2 x\ 
Zx-Z+l- — \y = 







a 2 
C 


±z , ^ 

— (xE- 

X 



E) 2 (4d) 

where dots denote differentiation with respect to the variable x. The particular rep- 
resentation of the Einstein-Maxwell system as given in (I4"al) - (l4d| ) may be easier to 
integrate in certain situation as demonstrated by Hansraj and Maharaj [17] and Ko- 
mathiraj and Maharaj 20| among others. To integrate the system (jlj) it is necessary 
to specify two of the variables y, Z, p, p or E. In our approach we choose Z and E 
on physical grounds. The remaining quantities are then obtained from the rest of the 
system 01]). 

In the integration procedure we make the choice 

Z= X -±^ (5) 
1 + x 

where k is a real constant. For the choice (jSJ) the line element (CQ) becomes 
ds i = -A 2 y 2 dt 2 + -J-±^—dx 2 + ^{d9 2 + sin 2 9d(f> 2 ) 

and we need to find the function y[x). The form (jSJ) ensures that the metric function 
behaves as 

e 2A = 1 + 0(r 2 ) 

near r = 0. In fact this is a sufficient condition for a static fluid sphere to be regular at 



the centre as pointed out by Maartens and Maharaj [23j. The choice (|5D ensures that 
the function is finite at the centre, remains regular and is well behaved in the stellar 
interior. This form contains particular solutions studied previously; for example when 
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k = — | we regain the neutron star model of Durgapal and Bannerji 21|. The choice (1ST) 
ensures that charged spheres found as exact solutions of the Einstein-Maxwell system 
will contain physically reasonable uncharged models when E = 0. 
On substituting dSJ) in (|4cl) we obtain 

E 2 (l+x) 2 ' 



4(1 + kx){l + x)y + 2{k - l)y + 



[1-k) 



Cx 



y = 



(6) 



which is the condition of pressure isotropy. The solution of the Einstein-Maxwell system 
(j3J, for the form §5§, depends on the integrability of fl§]). It is necessary to specify the 
electric field intensity E to integrate ([6]). A variety of choices for E is possible but only 
a few are physically reasonable which generate closed form solutions. We can reduce 
(JHJ) to simpler form if we let 

E?_ ax 

C (1 + x) 2 U 
where a is a constant. The form for E 2 in (J7|) vanishes at the centre of the star, and 
remains continuous and bounded in the interior of the star for a wide range of values 
of the parameter a. Thus the choice of E 2 is physically reasonable and useful in the 
gravitational analysis of charged spheres. When a = there is no charge, and we regain 
uncharged stellar solutions such as the neutron star model of Durgapal and Bannerji 



21] for the parameter value k = — |. Note that the same form for the electric field 



intensity (jTJ) was utilised by Hansraj and Maharaj 



17| | and Thirukkanesh and Maharaj 



18| . For this form of E 2 a comprehensive physical analysis was performed by Hansraj 



and Maharaj [17| for charged spheres when k = 0. They demonstrated that all criteria 
for physical acceptability were satisfied including the fact that the speed of sound is less 
that the speed of light. We expect that this form for the electric field intensity, with 
the generalised potential (jSJ), will provide a more general class of charged relativistic 
spheres with desirable physical features. 

Upon substituting the choice (j7j) in equation 



4(1 + kx)(l + x)y + 2(k - l)y 

which is the master equation for the system | 
introduced the transformation 



I we obtain 
[1-k- a)y 







(8) 



1 + x = KX, 

so that (jSJ) becomes 

AX(1-X) 



K 



k 



k 



Thirukkanesh and Maharaj 18] 



Y(X) = y{x) 



d 2 Y 
dX 2 



2§ + (K + a)Y 



0, 



a = — 



which is a special case of the hypergeometric equation. It is possible to integrate this 
hypergeometric equation. However it is important to note that the possible solutions 
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are restricted as k ^ and k ^ 1 because of the transformation used. We demonstrate 
in Section 3 that we can accommodate k = and k = 1 in a wider class of solutions. 



3 Master Equation 

It is convenient to introduce the new variable z = 1 + x in (jSJ) to yield 

Az(l - k + kzS^r + 2(k - 1)^$- + (l — k — a)y — (9) 

dz l dz 

where we have set y = y(z). Two categories of solution are possible for k — 1 and k ^ 1. 
Case I : A; = 1 

In this case ([2]) becomes the Euler-Cauchy equation with solution 



(l + x )(l+v^)/2 + C2 ( 1 + a; )(l-V^)/2 (10) 



where c\, c<i are constants. From (|4al) - (l4dj ), §J5§ and (11 01) we can show that the solution 
to the Einstein-Maxwell system, for our specific choices of Z and E, becomes 

e 2A = 1 (11a) 

e 2v = A 2 [ Cl (l + a;)( 1+ ^/ 2 + c 2 (l + i)' 1 -^)/ 2 ] 2 (lib) 
aCx 

P = -WTx? (UC) 



2C[ Cl (l + y^T^)(l + x)^ + c 2 (l- y^T^)] aCx 
P (1 + x)[ Cl (l + x)^ + c 2 ] 2(1 + xf { ' 

E - (lle) 

in terms of the variable x. Observe that this case is not regainable from Thirukkanesh 
and Maharaj [3] as k ^ 1 in their transformation. We do not pursue this case further 
as p < 0. 



Case II : k ^ 1 

As the point z = is a regular singular point of (Q, there exist two linearly indepen- 
dent solutions of the form of a power series with centre z = 0. These solutions can be 
generated using the method of Frobenius. Therefore we can write 

oo 

y = J2a t z l+ \ a ^0 (12) 

i=0 

where a, are the coefficients of the series and b is a constant. For a legitimate solution 
we need to determine the coefficients eij as well as the parameter b. On substituting 
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H in to (ED we obtain 



2a b(l-k)(2b-3)z b - 1 

oo 

J^{2(1 -k)(i + b + l)(2i + 26 - l)a m 



i=0 



+ [4fc(« + 6)(i + 6- 1) + (1 - A; - a)]^}^ 



b+i 







(13) 



in increasing powers of z. For equation (fT3"|) to hold for all powers of z in the interval 
of convergence we require 



2a 6(l-&)(26-3) = 



(14a) 



a-i+i 



4k{i + b)(i + b-l) + (1 - k - a) 
2(fc — l)(i + 6 + l)(2i + 26—1) 



a h i > (14b) 



Since ao 7^ and fc ^ 1, we have from f)14al) that 6 = or b — |. Equation (|14bl) is the 
basic difference equation governing the structure of the solution. It is possible to express 
the coefficient tij in terms of the leading coefficient ao by establishing a general structure 
for the coefficients by considering the leading terms. These coefficients generate the 
pattern 



H+l 



A Ak(p + 6) (p + b - 1) + (1 - k - a) 



p=0 



2(k - l)(p + 6 + l)(2p + 26 - l) 



(15) 



where we have utilised the conventional symbol H to denote multiplication. It is easy 
to establish that the result (TDj]) holds for all positive integers p using the principle of 
mathematical induction. 

Now it is possible to generate two linearly independent solutions to (Q with the 
help of (|12p and (fTBI) . For the parameter value 6 = we obtain the first solution 



Vi = ao 



i+EII 

i=0 p=0 



Akp(p - 1) + (1 - k - a) 
2(k- l)(p+ l)(2p- 1) 



J+i 



(16) 



For the parameter value 6 = | we obtain the second solution 



2/2 = a Q z 



3/2 



i=0 p=0 



fc(2p + 3)(2p+l) + (l-ib-a) 
(fc- l)(2p + 5)(2p + 2) 



(17) 



Thus we can write the general solution to the differential equation 
of the electric field given in (|7j), as 



for the choice 



V = c iVi + C22/2 



7 



where ci, c 2 are arbitrary constants, and y\ and §2 are given in (flEI) and (JTTj) respec- 
tively. In terms of the variable x we can write (|18[) as 



ci 



4fcjo(p - 1) + (1 — fc — a) 



1 + z 



i=0 p=0 

+ c 2 (l + a;) 3/2 



2(fc- l)(p+ l)(2p- 1) 

A fc(2p + 3)(2p+l) + (l-A;-a) m 
i + 2^11 (jfe-l)(2p + 5)(2p + 2) 1 J 



i=0 p=0 



ciyi(x) + c 2 t/ 2 (x) 



(19) 



where we have set c\ = ciclq and c 2 = c 2 ao for simplicity. The solution to the Einstein- 
Maxwell system fl4ap - (l4dl) . for our specific choices of Z and E, can now be written as 



e 2A 

e 2v 

P_ 

C 

P_ 

c 
&_ 
c 



x 



1 + kx 
A 2 y 2 

(1 - k)(3 + x) 
(1 + x) 2 



ax 



2(1 + x) 



i (l + kx)_ y | (fc-1) | 



a x 



[1 + x) y 

ax 



1 + x) 2(1 + x) 2 



1 + x) 



(20a) 
(20b) 
(20c) 

(20d) 

(20e) 



in terms of the variable x. The form of the exact solution (120a[) - fl20ej) has a similar 
structure to the general solution of Thirukkanesh and Maharaj 181 ]; however it is 
important to realise that our solution is a new result because the series in (fl~9l) is 
different. In addition note that k = is allowed in fl20al) - fl20el) unlike the result of 
Thirukkanesh and Maharaj 3]; our result can be interpreted as a generalisation. 



4 Elementary solutions 

The general solution ffl9|) is given in the form of a series which define special functions. 
For particular values of the parameters involved it is possible for the general solution 
to be written in terms of elementary functions which is a more desirable form for the 
physical description of a charged relativistic star. We find two linearly independent 
solutions, in terms of elementary functions, in this section for the differential equation 
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4.1 First solution 



On substituting 6 = and setting 1 — k — a = —Akn(n — 1) in (|14bl) we obtain 

4k (n — i)(n + i — 1) , s 

Oi+i = ^tt^ i > 21 

+ 1-k (2i + 2)(2i-l) V ; 

It is easy to see from ( 12TT) that a n+ i = 0. Clearly the subsequent coefficients a n+ 2, a ra +3, 



a n+ 4, . . . vanish. Equation (Izll) has the solution 

4fc \ i n(n-l)(2i-l)(n + i-2)! 



a,- = - 



1 - k 



a>o, < i < n 



(2i)\(n-i)\ 

Then from f|T2|) (when 6 = 0) and fl22|) . we obtain 

4fc Vn(n- 1)(2«- l)(n + i- 2)! 



(22) 



yi = «o 



k 



(2i)\(n-i)\ 



(23) 



for 1 — k — a = —4kn(n — 1). 

On substituting b = | and setting 1 — /c — a = —k(2n + 3)(2n + 1) in (114b|) we 
obtain 

4A; (n-z)(ra + z + 2) 



-Oj, 2 > 



(24) 



1-k (2i + 5)(2i + 2) 
We observe from ( 1241) that a n+ i = 0. Consequently the remaining coefficients a n+ 2, a n+ s, 
a n+ 4, . . . vanish. Equation (|24l) may be solved to yield 



4k V (2z + 2)(n + z + l)! 

(n + l)(2i + 3)!(n-i)! 



Now from (fT2|l (when 6 = |) and (j25|) . we obtain 



ao, < z < n 



(25) 



3/2 



(2i + 2)(n + z + l)! 
fcy (n+l)(2i + 3)!(n-i)!' 



(26) 



for 1 - Jfc - a = -k(2n + 3)(2n + 1). 

The polynomial and algebraic functions (|23l) and (j26|) comprise the first solution to 
the differential equation ([9]) for appropriate values of k and a. 



4.2 Second solution 

We can find the second solution of (Q using the method of reduction of order in 
principle. However this proves to be difficult in practice because of the complicated 
form of the first solution given in (123]) and (126]) . We utilise a transformation to first 
simplify (jHJ) before seeking the second solution. We take the second solution of (jHJ) to 
be of the form 

y = (1 - k + kz) 1/2 u{z) (27) 
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when u(z) is an arbitrary polynomial. Special cases of ([27)1 are known to solve (Q 
which motivates the algebraic form for y as a generic second solution to the differential 
equation On substituting y in <$f) we obtain 

4z(l -k + kz)j^ + [Akz + 2(k - 1)]— + (1 - 2k - a)u = (28) 

which is linear differential equation for u(z). 

As in Section 4.1 it is possible to find solutions in terms of polynomials, and product 
of polynomials with algebraic functions for u(z) for certain values of the parameters k 
and a. As the point z = is a regular singular point of (128]) . there exist two linearly 
independent solutions of the form of the power series with centre z = 0. Thus we 
assume 

£>* <+d (29) 



oo 

U 

i=0 

where the constants c,- are the coefficients of the series and d is the constant. Substi- 



tuting (T5HD in (J2SJ) we obtain 

2c d(l- k)(2d-3)z d - 1 

oo 

+ ^{2(l-fc)(i + d+l)(2« + 2d-l)ci+i 

i=0 

+ [4fc(i + df + (1 - 2k - a)}a}z i+d = (30) 

For equation (130|) to hold true for all z we require that 

2c d(l- k)(2d-3) = (31a) 

Ak{i + df + {l-2k-a) 
Ct+1 = 2{k-m + d + l){2r + 2d-l) C - Z -° (31b) 

Since Co 7^ and fc 7^ 1 we have from (I31al) that d — or d — |. Equation (I31bl) is the 
linear recurrence relation governing the structure of the solution. 

On substituting d = and setting 1 — k — a = —k(2n + l)(2n + 3) in (131bj) we 
obtain 

Q4-1 = r ; w ; Cj, « > 32 

+ \l-fc/ (2i + 2)(2i-l) v ; 

From K32J we have that c n +2 — and subsequent coefficients c n _|_3, c n +4, c n +5, . . . 

vanish. Then ( l32l) has the solution 

/ 4A; V (n + l)(2i- l)(n + i)! , s 

C ' = - (t^i) ( g l (n-« + l)l V ° S 12 " + 1 (33) 



Then from (129]) (when d = 0) and (133]) we obtain 



V2 = C 



ra+1 



i=l 



4k V (n+l)(2«-l)(n + i)! 



(2i)!(re-i + l)! 



-2 



(34) 
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for 1 - k - a = -k(2n + l)(2n + 3). 



On substituting d — § and setting 1 — k — a = —4kn(n — 1) in (131bl) we obtain 



Ci+l 



4A; \ (n + i + l)(n -i-2) 



Ci, i > 



(35) 



(2i + 2)(2z + 5) 

From (I3"5"j) we have that c n __i = and subsequent coefficients c n , c n+ i, c n+2 , . . . vanish. 
The equation (|35|) has the solution 

4k V {2i + 2){n + i)\ 



l-kj n{n- l)(2i + 3)!(n-i-2)! 



c , < i < n — 2 



From (1291) (when <i = |) and (J2SJ) we obtain 



n-2 



'^1 

i=l v 



4fc 



{2i + 2){n + i)\ 



(36) 



(37) 



A;/ n(n- l)(2i + 3)!(n -i-2)! 

for 1 — k — a = —4kn(n — 1). 

The algebraic solutions (l34j) and (1371) comprise the second solution of the differential 
equation Q) which are clearly independent from (I26I) and (1231) . 



4.3 Elementary functions 

We have obtained one class of polynomial solution ( l23l) and three classes of solutions 
( 1261) . ( 1341) . and ( 1371) in terms of product of polynomials and algebraic functions. The 
polynomial solution (123]) and the product of polynomials with an algebraic function 
(I26I) generate the first solution. The second linearly independent solution is given by 
( |34l) and ( !37l) which are products of polynomials and algebraic functions. By collecting 
these results we can express the general solution to (GO) in two categories. We express 
the general solution in terms of the independent variable x. From (1231) and (1371) we 
write the first category as 

4k y n(n-l)(2i-l)(n + i-2)! < 
-k) (2i)\(n-i)\ K ' 

n-2 



y 



A 



i=l v 



+ £(l + fcr) 1/2 (l+:r) 3/2 



3 Mi 

i=l v 



4fc 



(2z + 2)(n + z)! 



-k) n(n- l)(2i + 3)!(n-z-2)! 



(1 + x) 1 



(38) 



for 1 — A; — a; = —4kn{n — 1). From (1261) and ( 1341) . the second category of solution has 
the form 

4k 



y = A(l + xf 2 
+ B(l + kx) 1/2 



n.+l 



(2z + 2)(n + z + l)! 
-kj (n + l)(2i + 3)\(n - i)\ 



i=i v 



1/,- 



-A; 



(n + l)(2i - l)(n + i)\ 
(2i)!(ra-i + 1)! 



1 + x 



1 + x) 



(39) 



11 



for 1 — k — a = — k(2n + l)(2n + 3). In the above A and B are arbitrary constants. Con- 
sequently we have demonstrated that elementary functions can be extracted from the 
general series f)19p by restricting the parameter values a and k. The general solutions 
(I38I) and (|39|) have a very simple form. It is important to observe that the Einstein- 
Maxwell solutions fl38l) and fl39l) apply to both charged and uncharged relativistic stars. 
We regain neutral exact solutions, which may be possibly new, by setting a = 0. 



5 Known cases 

We may generate individual models for charged and uncharged stars found previously 
from our general class of solutions. These can be explicitly regained from the general 
series solution (T0?|) or the elementary functions (1551) and (1591) . We demonstrate that 
this is possible in the following cases: 



Case I : Thirukkanesh and Maharaj charged stars 

If we set K = ^=i, a — j, C\ = di and c 2 = K~ 3 / 2 d 2 then (|TS1) can be written as 



di 



+ d 2 



1 + x 
K 



oo i 

En 

3/2 



Ap(p - 1) - (K + a) fl + x 
2(p+l)(2p-l) V K 



i+l 



oo i 
i=0 p=0 



(2p + 3)(2p+ 1) - (K + a) fl + x 



(2p + 5)(2p + 2) 



i+l 



(40) 



The Einstein-Maxwell solution (I4UI) corresponds to the charged model of Thirukkanesh 
and Maharaj [3]. However note that k and A; 7^ 1 in fj4l7]) . In our wider class of 
solutions f|T9l) it is permitted that k = 0; the exact solution ffTUj) corresponds to the case 
k — 1. The Thirukkanesh- Maharaj charged stars contain neutron star models found 
previously including the Durgapal and Bannerji model, and they are consequently 
physically reasonable. 



Case II : Hansraj and Maharaj charged stars 
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If we set k = 



then it is possible after some manipulation, to write (fT9|) in the form 



Cl 



+ c 2 {l + x 



i=0 
\3/2 



^«S±!| (1 -.,, 1+If 



i=0 



(2i + 5)! 



Cl 



E 



i=0 



1 



1 



+ 



3c 2 



a 



)3/2 



1 oo 

3-D 



i=0 



(2^ + 2)! (2^ + 1)! 
[v/(l- a )(l + x)] 3 

1 1 

(2i + 5)! ~ (2i + 4)! 



h/(l-a)(l + s)] 2 ( i+1 > 



(41) 



It is interesting to observe that the last equation above can be expressed in terms of 
trigonometric functions. Then it is easy to show that (14 ip can be written in the form 

V = C a [COS - a)(l+x) + y/(l - Ct)(l + x) SU1 y/ (1 - Ct)(l + x)] 



+ 



3c 2 



I sin a/(1 - a)(l + x) - a/(1 - a)(l + a;) cos a/(1 - a)(l + x)] 



(42) 



is 



(1 - a) 3 / 21 

= [di + d 2 y/(l - a)(l + x)] sin ^(1 - a)(l + x) 
+ [d 2 - di v 7 ! 1 - «)(! + z)] cos V(l - a)(l + x) 

where we have set d\ = ^^3/2 an d d 2 = C\. The Einstein-Maxwell solution 
the same as the charged model of Hansraj and Maharaj 17]. Note that our solution 
(1421) corrects a minor misprint in the result of 17J . The Hansraj-Maharaj charged stars 
were comprehensively studied in 17J; it was demonstrated that their model produced 
a charged relativistic sphere that satisfied all physical criteria. In particular the speed 
of sound was less that the speed of light and causality is maintained. 



Case III : Finch and Skea neutron stars 

If we set k = and a = 0, and follows the procedure outlined for Case II, then f[T$l) 
becomes 

y = [dx + d 2 y/(l + x)] sin y/(l + x) + [d 2 - d 1 y/(l + x)] cos y/(l + x) (43) 

where we have set d\ = 3c 2 and d 2 = C\. Alternatively we can obtain the result (j43p 
directly from fT42"]) w ith a = 0. The exact solution (I43p is the neutron star model of 
Finch and Skea |22']. The Finch-Skea model satisfies all the physical conditions for 
an isolated spherically symmetric stellar source, and consequently has been utilised by 
many researchers to model neutron stars. 
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Case IV : Durgapal and Bannerji neutron stars 

If we set a = and k = — ~(n = 0), then (15^1) becomes 



y = c(l + xf /2 + d(2 



x 



,1/2 



(5 + 2x) 



(44) 



where we have set c = A and d = The exact solution (1441) was first found by 

Durgapal and Bannerji 21]. The Durgapal-Bannerji solution has been widely applied 
as a relativistic model for neutral stars with superdense matter. 



6 Discussion 

We have found new solutions (1201) to the Einstein-Maxwell system (T4J) by utilising the 
method of Frobenius for an infinite series; a form for one of the gravitational potentials 
was assumed and the electric field intensity was specified. These solutions are given in 
terms of special functions. For particular values of the parameters involved it is possible 
to write the solution in terms of elementary functions: polynomials and products of 
polynomials and algebraic functions. The electromagnetic field may vanish in the 
general series solutions and we can regain uncharged solutions. Thus our approach 
has the advantage of necessarily containing a neutral stellar solution. The charged 



model of Thirukkanesh and Maharaj 18j, the charged spheres of Hansraj and Maharaj 



171 ]. the uncharged neutron star model of Finch and Skea [22], and the uncharged 
superdense star of Durgapal and Bannerji {2l| are contained as special cases in our 
class of solutions to the Einstein-Maxwell system. The simple form of the solutions 
found facilitates the analysis of the physical features of a charged sphere. 

We are now in a position to investigate graphically the behaviour of the matter 
variables in the stellar interior for particular choices of the parameter values in (120]) . 
We have generated figure 1 with the assistance of the software package Mathematica. 
For simplicity we make the choices A — B — C — 1, k— j and a — |, over the 
interval < r < 1, to generate the relevant plots. In figure 1, plot A represents the 
energy density p, plot B represents the pressure p, and plot C represents the electric 
field intensity E 2 . It can easily be seen that these matter variables remain regular in 
the interior of the star for < r < 1. We observe that the energy density and the 
pressure are positive and monotonically decreasing functions in the interior of the star. 
The electric field intensity is positive and monotonically increasing for the interval of 
the plot. Thus the quantities p, p and E 2 are continuous, regular and well behaved in 
the interior of the star. Consequently the solutions found in this paper are likely to be 
useful in the description of charged relativistic fluid spheres. 
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Figure 1: Plots of the energy density p, pressure p and electric field E . 
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